Collective Josephson vortex dynamics in a finite number of intrinsic Josephson 
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We report the experimental confirmation of the collective transverse plasma modes excited by 
the Josephson vortex lattice in stacks of intrinsic Josephson junctions in BiaS^CaC^Os+a single 
crystals. The excitation was confirmed by analyzing the temperature (T) and magnetic field (H) de- 
pendencies of the multiple sub-branches in the Josephson-vortex-flow region of the current-voltage 
characteristics of the system. In the near-static Josephson vortex state for a low tunneling bias 
current, pronounced magnetoresistance oscillations were observed, which represented a triangular- 
lattice vortex configuration along the c axis. In the dynamic vortex state in a sufficiently high mag- 
netic field and for a high bias current, splitting of a single Josephson vortex- flow branch into multiple 
sub-branches was observed. Detailed examination of the sub-branches for varying H field reveals 
that sub-branches represent the different modes of the Josephson-vortex lattice along the c axis, with 
varied configuration from a triangular to a rectangular lattices. These multiple sub-branches merge 
to a single curve at a characteristic temperature, above which no dynamical structural transitions 
of the Josephson vortex lattice is expected. 

PACS numbers: 74.72.Hs, 74.50.+r, 74.78.Fk, 85.25.Cp 



I. INTRODUCTION 



pressed as 



-a 

G 

O 

o 



> 

m 
m 
O 

o 

-I— > 
c3 



i 

C 

O 

o 



X 



It has been known that stacked intrinsic Joseph- 
son junctions (IJJs) form in naturally grown 
Bi 2 Sr 2 CaCu 2 8+I (Bi-2212) single crystals,! The 
collective Josephson plasma oscillation, manifested 
by the electromagnetic resonant absorption in such 
Josephson-coupled layered superconductors, has pro- 
vided the key to understanding the superconducting 
state and the vortex phases^ forming in the materials. 
Two kinds of collective Josephson plasma modes, lon- 
gitudinal and transverse, exist in such a system of an 
infinite number of staked IJJs£ The longitudinal plasma 
mode oscillation, excited in an externally applied o 
axis-oscillating electric field, propagates along the c axis 
of a stack of IJJs with the inter-layer phase difference 
of a junction being unform in a planar direction. The 
resonant absorption of externally applied microwaves 
by the excitation of the longitudinal modes enables one 
to directly obtain the Josephson plasma frequency, u> p . 
This, in turn, provides information on the interlayer 
Josephson coupling strength and the c-axis superfluid 
density in Bi-2212. 4 By contrast, the transverse plasma 
modes can be excited by the oscillating magnetic fields in 
the ab plane. The resulting induced current along the c 
axis generates the oscillating electric field that resonates 
with the transverse Josephson plasma oscillation. 

On the other hand, a system with a finite number 
(N) of stacked IJJs, described by linearized coupled sine- 
Gordon differential equations with the inductive inter- 
junction coupling, exhibits collective transverse plasma 
(CTP) eigen modes, with the eigen frequencies lying be- 
tween those of the longitudinal and the transverse modes 
found for a system of an infinite number of stacked 
junctions.^ The corresponding dispersion relation is ex- 
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Here, A a & (A c ) is the c-axis (afr-plane) London penetra- 
tion depth, N the number of the junctions in a stack, 
s=0.3 nm (_D=1.2 nm) the thickness of a superconducting 
(insulating) layer, fc(=27rm/L; L=the length of a stack; 
m=the number of the vortices present in a junction) the 
wave number of the plasma oscillation^ The mode index 
n runs from 1 to N. 

The CTP modes can be excited by the moving Joseph- 
son vortex lattice forming in high magnetic fields, which 
are alternative solutions of the coupled sine-Gordon 
equations^ The temporal oscillation of the inter-layer 
phase difference due to the moving Josephson vortex 
lattice, with the frequency matching with that of any 
CTP modes, excites a resonant plasma oscillation. The 
Josephson vortex lattice also transforms its spatial con- 
figuration along the c axis in accordance with the c-axis 
standing-wave modes of the plasma oscillation. The reso- 
nance of Josephson vortex lattice to the CTP modes is re- 
vealed in the form of the multiple Josephson vortex-flow 
branches in the tunneling current- voltage (I-V) charac- 
teristics of a stack of IJJsAS 

Experimental observation of the multiple CTP modes 
in a stack of a finite number of IJJs has often been at- 
tempted in a mesa or an S-shaped stack of IJJs fabricated 
on the surface of a large stack of IJJs, which is called a 
pedestal or a basal layeri 10 : 11 Since each superconduct- 
ing bilayer in a stack of IJJs is much thinner than the 
c-axis London penetration depth A a b, Josephson vortices 
in a usual mesa or an S-shaped-stack structure are in a 
strong inductive coupling to those in the basal layer that 
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FIG. 1: The phase diagram of the Josephson vortex configu- 
ration for different ranges of transverse magnetic field and the 
length of the intrinsic- Josephson-junction stack. A schematic 
figure in the each region represents the supercurrent distribu- 
tions in layered structures, where the centers of the Josephson 
vortices are indicated by black dots. 




contains a large number of IJJs by itself, although mea- 
surements are intended to be made on a finite number 
of stacks in the structure of interest. Thus, the stack of 
IJJs with a coupling to the basal layers corresponds to 
a system containing nearly an infinite number of junc- 
tions, which does not support the CTP mode formation 
along the c axis. This situation significantly affects the 
dynamics of the Josephson vortices in the mesa or the 
S-shaped stack itself. This may explain the usual fail- 
ure of observing the multiple resonating modes in those 
structures. 



In this paper, we present the static and dynamic states 
of Josephson vortices in stacks of IJJs, each sandwiched 
between two Au electrodes without the basal layer, thus 
consisting of a genuinely finite number of IJJs. In the 
near-static state in a low tunneling bias current on a 
single Josephson vortex-flow branch, pronounced mag- 
netoresistance oscillations were observed with h/Ae mag- 
netic flux period. This half-flux-quantum periodic resis- 
tance oscillations indicate that the Josephson vortex lat- 
tice in the near-static state constitutes a triangular lattice 
along the c axis. In the dynamic state in a high tunnel- 
ing bias current, on the other hand, a single Josephson 
vortex-flow branch splits into multiple sub-branches. De- 
tailed examination of the sub-branches for varying mag- 
netic fields and temperatures revealed that the observed 
sub-branches were caused by the resonance between the 
Josephson vortex lattices and the CTP modes. 



II. BASIC CONCEPTS 

The Josephson vortex configuration in stacked IJJs of 
Bi-2212 single crystals is determined by the magnetic flux 
density in the junctions and the junction sizeil Joseph- 
son vortices are introduced to the insulating layer of a 
junction when its length L is longer than the Joseph- 
son penetration depth, Xj {=y / &o/2Trn j c (D + 2A 2 b /s)], 
where $o(=^-/2e) is a flux quantum. One can identify 
various characteristic magnetic fields in relation with the 
dynamics of Josephson vortices [see Fig. 1]. For a stack 
of junctions of length L Josephson vortices start enter- 
ing the junctions above a characteristic fieldji&ii 
[=($o \ c /L 2 \ a b)ln(X a b /D)] . This field can be higher than 
the nominal value of one flux quantum threading a junc- 
tion of length L, that is (=$ /(s + D)L), especially 

for a junction comparable to Aj. Even above H]) the 
tunneling screening current prevents an external mag- 
netic field from entering the junctions until it reaches 
hJl- Thus, for a junction with its length in the range of 
Xj < L < L m [=Ajln(A a b/(s + D)], one can observe the 
usual oscillatory Fraunhofer modulation of the tunneling 
critical current as a function of an external magnetic field 
of H < h]?!, although the junction size is larger than Xj 
(Region 0)J& 

The Josephson vortices entering the junctions for H > 
flj.j form a regular lattice as the magnetic field reaches 

the valued (=1.4$ /27r(s+£»)Aj). In this field range 
(Region 1), since the field-induced vortices exist in arbi- 
trary configurations in a given magnetic field, disordered 
multiple vortex-flow subbranches are generated in the I- V 
characteristics. In this region, the number of vortices in 
a junction and the propagation velocities of different vor- 
tex configurations are not correlated with each otherii£ 
Thus, in this field range, the vortex-flow subbranches in 
the I- V characteristics are not related to vortex mode 
velocities&i& or geometrical resonances^ 

In an external field above Hj} 2 (Region 2), Josephson 
vortices start forming a triangular lattice. Recent ob- 
servation of the oscillatory vortex-flow resistance aris- 
ing from the interaction between the vortex lattice and 
the crystal boundary potential clearly confirmed the ex- 
istence of the triangular vortex lattic e) 20 ! 21 above h]} 2 . 

On the other hand, the resonance between the Joseph- 
son vortex lattice and CTP modes is well known in 
this vortex dynamics system^ The collective resonance 
revealed as voltage peaks (in the case of the voltage 
bias) or voltage jumps (in the case of the current bias) 
has been searched extensively both theoretically and 
experimentallyi£ii2i22 The collective resonance has been 
observed in the dense-vortex state corresponding to the 
field range of H > fzjg [=$o/2A, 7 (.s + £>)], where the 
shortest inter-vortex spacing becomes comparable to the 
diameter of a Josephson vortex 2Aj (Region 3). Thus, 
Josephson vortices in Regions 2 and 3, although they 
form triangular lattices in both regions in a static state 
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without a tunneling bias current (Fig. 1), behave differ- 
ently in the vortex-flow state in the presence of a finite 
tunneling bias current. In the dynamic state, Joseph- 
son vortices in Region 2 remains in a triangular lattice, 
while those in Region 3 transform their lattice configura- 
tion to fit the plasma propagation modes in the stacked 
junctions^ The non-Josephson-like emission from the 
Josephson-vortex motion^ and the Shapiro resonance 
steps in the Josephson vortex-flow stated have been ob- 
served in this dense-vortex state. These results confirm 
the existence of the coherent motion of the Josephson 
vortex lattice through the entire thickness of stacked (in- 
trinsic) Josephson junctions. 

For the highest characteristic field [=$ /Aj(s + Z?) 
(~4 T for Bi-2212)], the shortest inter- vortex spacing be- 
comes comparable to the Josephson penetration depth, 26 
A,/. In this region (Region 4), the Josephson current 
along the length of a junction distributes almost sinu- 
soidally so that the collective resonance behavior between 
the vortex motion and the plasma modes is expected to 
be stronger than in a field lower than H^. 



III. EXPERIMENT 

As-grown slightly overdoped Bi-2212 single crystals 
were prepared by the conventional self-flux method. We 
fabricated, using the double-side cleaving technique, 27 
three samples of IJJs sandwiched between two Au elec- 
trodes at the top and the bottom of a stack of IJJs with- 
out the basal part [Fig. 2(b)]. This structure, which is in 
contrast to the usual mesa or the S-shaped-stack struc- 
ture fabricated on the surface of a single crystal with a 
large basal layer, is more coincident with the model sys- 
tems usually adopted in the theoretical analysis. 18 The 
geometry with the basal layer eliminated by the double- 
side cleaving enables one to investigate the Josephson 
vortex dynamics in coupled IJJs without the interference 
from the vortex motion in the basal layer. Furthermore, 
the top and the bottom Au-film leads in a stack help 



the formation of the standing plasma oscillation modes 
along the c axis, while forming a transmission guide of 
the plasma oscillations along the length of the stack. 

The samples containing sandwiched stacks were pre- 
pared in the following way^£ A single crystal was first 
glued on a sapphire substrate using negative photoresist 
and was cleaved until an optically clean and flat sur- 
face was obtained. Then, a 100-nm-thick Au film was 
thermally deposited on top of the crystal to protect the 
surface of the crystal and to make a good ohmic con- 
tact. A mesa structure was then prepared by photolitho- 
graphic patterning and Ar-ion-beam etching. The surface 
of the patterned mesa was fixed to another sapphire sub- 
strate using negative photoresist and the basal layer of 
the mesa was subsequently cleaved away. This process, 
called the double-side-cleaving technique, was developed 
by Wang et al2L A 100-nm-thick Au film was again de- 
posited on this freshly cleaved crystal surface, leaving a 
stack of IJJs sandwiched between two Au electrodes. A 
few-/um-long portions on both sides of the sandwiched 
stack were etched away subsequently [see Fig. 2(a)] to 
get the bottom Au electrode exposed for c-axis transport 
measurements as schematically depicted in Fig. 2(b). Fi- 
nally, a 300-nm-thick Au-extension pad was attached on 
each of the Au electrode. The lateral sizes of three sam- 
ples were 13.5 x 1.4 fim 2 [SP1], 15 x 1.4 ^m 2 [SP2] and 
16 x 1.5 fim 2 [SP3]. 

The magnetic field was aligned in parallel with the 
plane of junctions within the resolution of 0.01 degree to 
minimize the pinning of Josephson vortices by the pan- 
cake vortices that can form in CuC>2 bilayers due to the 
vertical component of a misaligned field. The alignment 
was done, in a field of 4 T and at a temperature around 
60 K, by tuning to the angle that gives the maximum 
Josephson-vortex-flow resistance or the minimum pining 
of the Josephson vortices by the pancake vortices gener- 
ated by the angle misalignment^ To minimize the exter- 
nal noise, transport measurements were performed with 
a low-pass filter connected to each electrode located at 
room temperature. 




FIG. 2: (color online) (a) An optical micrograph of a sample 
during the fabrication process. The scale bar represents the 
length of 10 fim. (b) Schematics of sample and measurement 
configurations. 



TABLE I: Sample parameters. TV is the total number of in- 
trinsic Josephson junctions in a stack, Aj the Josephson pen- 
etration depth, jyj.3=$o/2Aj(s + D) (s and D the thickness of 
the superconducting and the insulating layers, respectively), 
Cm the maximum propagation velocity of vortices in Region 
3, and ci the velocity of the fastest collective plasma mode 
estimated by the RCSJ model. 
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FIG. 3: (color online) (a) The Josephson vortex-flow resis- 
tance, R as a function of an external magnetic field H for 
SP1 at T—55 K in bias currents from 1 /iA to 10 /iA at inter- 
vals of 1 /iA. (b) Temperature dependence of the Josephson- 
vortex-flow sub-branches and quasiparticle branches of SP1 at 
if =3. 9 T. Inset of (a): a single Josephson vortex-flow branch 
and quasiparticle tunneling branches in H—1.2 T at T—55 K. 
Inset of (b): the temperature dependence of R(H) in a bias 
current of 3 /iA. Curves in the main panel and the inset of 
(b) are shifted vertically for clarity. 



IV. RESULTS AND DISCUSSION 

A. Static and dynamic state 

Fig. 3(a) shows the oscillation of Josephson-vortex- 
flow resistance of SP1 [7V=35] as a function of external 
magnetic fields at T=55 K in bias currents from 1 /iA 
to 10 /iA at intervals of 1 /iA. The field period of the 
oscillation, H p =hl0 G is one half of H (=$ /£(s + D)) 
through the junction area of length L=13.5 /im and thick- 
ness (s + D)=1.5 nm. The oscillation of the Josephson- 
vortex-flow resistance is known to be caused by the in- 
teraction between the triangular Josephson vortex lattice 
along the c axis and the vortex-flow boundary poten- 
tial at the junction edges ^LSLSSi The observation of the 
Josephson-vortex-flow resistance oscillation in our sand- 



wiched stack without the basal layers indicates that, as 
the S-shaped-stack geometry, our sandwiched geometry 
also supports the Josephson-vortex triangular lattice in 
the steady state [see the inset of Fig. 4(a)]. The oscil- 
lation in Fig. 3(a) is reduced by increasing bias current 
as the bias current tilts the boundary potential. 3 - The 
oscillatory-magnetoresistance behavior is also suppressed 
completely at T ^80 K near T c (=92 K) as in the inset of 
Fig. 3(b), which is related to melting of the near-static 
Josephson vortex latticed 

The inset of Fig. 3(a) shows the single Josephson- 
vortex-flow branch (the low-bias region) and quasiparti- 
cle branches (the high-bias region) in H—1.2 T at T=55 
K. The field, belonging to Region 2 in Fig. 1, corresponds 
to the position indicated by the arrow in the main panel. 
In this field region, the single-curve Josephson-vortex- 
flow branch caused by the triangular lattice persists even 
down to 4.2 K (not shown). In Fig. 3(b), /-^character- 
istics in a field of Region 3 at T=25 K also show a single 
Josephson vortex-flow branch corresponding to the trian- 
gular lattice. Table I shows the characteristic field value 
of H^ of each sample. 

Contrary to the field range of Region 2, however, the 
single Josephson-vortex-flow branch in Region 3 evolves 
into the multiple Josephson-vortex-flow sub-branches 
with decreasing temperature. Fig. 3(b) shows the 
temperature dependence of the Josephson-vortex-flow 
branches, where the multiple sub-branches start appear- 
ing at T^20 K and become clearer with decreasing tem- 
perature. At T=4.2 K, a single branch representing the 
triangular lattice below the current bias of /b~2 /iA splits 
into multiple Josephson-vortex-flow sub-branches for the 
higher current bias. The Josephson vortex flow, then, 
transits to the McCumber quasiparticle tunneling for /{, 
above 8 /iA. The multiple Josephson-vortex-flow sub- 
branches rapidly suppress with increasing temperature 
above 4.2 K and disappear completely as the tempera- 
ture reaches 25 K. On the contrary, the quasiparticle- 
tunneling branches are robust to the temperature vari- 
ation and persist up to T=25 K. This contrasting tem- 
perature dependencies of two different regions strongly 
imply that the multiple branches in the two regions are 
of different physical origins. 

We also observed the separation of the two charac- 
teristic multiple-branch regions in SP2 and SP3. Figs. 
4(a) and 4(b) show the magnetic field dependencies of 
/- V characteristics of two sandwiched stacks, SP2 and 
SP3, respectively, at 4.2 K. The number of zero-field 
quasiparticle-tunneling branches of SP2 and SP3 were 22 
and 45, respectively (not shown), which coincided with 
the total number of IJJs in the given stacks. The differen- 
tiation of the multiple branches by the collective Joseph- 
son vortex flow (the low-bias region) from the ones by 
the quasiparticle tunneling (the high-bias region) is based 
on the fact that the delimiting current, represented by 
the current values marked by the dotted vertical lines, 
stemmed from the zero- field Josephson critical current. 
Nonetheless, a finite resistance is present in the region 
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FIG. 4: (color online) The magnetic field dependencies of 
the Josephson vortex-flow sub-branches and the McCumber 
branches of (a) SP2 and (b) SP3 at T=4.2 K. The I-V char- 
acteristics in all the figures are shifted vertically for clarity. 
A voltage V m denoted by a vertical dotted line in each set of 
curves demarcates the Josephson vortex flow region (JVFR) 
from the McCumber quasiparticle-tunneling region (QPR). 
Inset of (a): (i) the triangular lattice and (ii) the rectangu- 
lar lattice configurations of the Josephson vortex lattice along 
the c axis. The contact resistance was subtracted numerically 
in the main panels of (a) and (b). 



below the critical bias point because of the Josephson 
vortex flow, although the region itself represents the dis- 
sipationless pair tunneling state. 

The observed sub-branches for SP2 and SP3 have com- 
mon features: (i) a single Josephson-vortex-fiow branch 
splits into multiple sub-branches for Regions 3 and 4, 
(ii) the number of the quasiparticle sub-branches is ~18 
and ^42, respectively, which is close to the number of 
the junctions in each sample, and (hi) the sub-branches 
in the Josephson vortex-flow region become clearer and 
wider with increasing the transverse magnetic field (in 
contrast to the case of quasiparticle sub-branches, which 
keep shrinking with field instead). These three common 
properties of sub-branches are consistent with the pre- 
diction of the inductive coupling theory for the electro- 



dynamics of stacked Josephson junctions. 

Since the speed of vortices, c m , is limited by that of the 
maximum electromagnetic wave in a junction, Josephson- 
vortex-fiow state changes to quasiparticle tunneling state 
near the maximum cut-off velocity or equivalently by the 
limiting voltage V m in the outermost sub-branch. Using 
the relation^ of V m =NH(s + D)c m , we could get c m as 
~3.9xl0 5 m/s and ^4.2 xlO 5 m/s for the correspond- 
ing delimiting bias point V m in SP2 and SP3, respec- 
tively. These values are similar to those obtained from 
a Josephson-vortex-fiow branch in mesa structures^ In 
the case of the usual mesa or the S-shaped-stack struc- 
ture with a basal layer, only a single Josephson-vortex- 
fiow branch used to be observed below the critical bias 
point, regardless of the temperature and the magnetic 
field intensity. The velocity c m corresponding to the de- 
limiting value on this single branch has been regarded as 
the slowest mode velocity, n=N. As illustrated below, 
however, in our case without the basal layer c m , corre- 
sponding to V m in the outermost sub-branch, represents 
the propagation velocity of the n=l plasma mode. 

In the resistively and capacitively shunted junction 
(RCSJ) models the junction capacitance is obtained by 
the relation, Cj—f3 C 2 ^ a R:1 , where (3 C [—(4I c /nI r ) 2 ] is 
the McCumber parameter, R n is the normal-state junc- 
tion resistance and I r is the return current. In SP2, 
Cj is obtained to be about 47 pF, with J c (=0.2 mA), 
I r (=4.8 fiA) and i?„(=10 tt). The Swihart velocity 
co=y/sA/2fiQCjX^ b , estimated from the junction capac- 
itance, is 3.6 x 10 4 m/s, where A o b=200 nm and the 
junction area of A=21 /im 2 . The resulting fastest mode 
velocity Cx=3.7 x 10 5 m/s is in remarkable agreement 
with the observed value of c rn in SP2. Experimental val- 
ues of c m in Table I are consistent with the calculated 
ones of Ci for all three samples. This indicates that the 
outermost sub-branch is the fastest velocity, n=l, mode. 

In general, Josephson-vortex resonance branches by 
the inter-layer inductive coupling have been predicted 
to generate current steps near the resonance voltages 
on the background of a single resistive branch^ In the 
present study, however, Jospehson-vortex resonance to 
the CTP modes constitutes multiple sub-branches over 
almost the entire bias range. A recent theory for an 
interlayer inductive coupling hybridized by a capacitive 
coupling predicts that the resistive sub-branches corre- 
sponding to the collective resonance modes show linear 
/- V characteristics in the low-bias region with the var- 
ied slopes from one another^ More spread in the slopes 
are predicted for a higher capacitive coupling. Fig. 5(a) 
shows the Josephson-vortex-fiow sub-branches for SP2 at 
H=5 T. As predicted by the theory, all the sub-branches 
in the Josephson vortex-flow regime exhibit linearities 
with varied slopes in the low bias range. The extrap- 
olation of the linear low-bias region converges to a sin- 
gle point on the current axis, although the inductive- 
capacitive hybrid coupling theory predicts they should 
converge to the origin. The discrepancy indicates a fi- 
nite pinning of Josephson vortex motion due to the pres- 
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FIG. 5: (a) Josephson vortex-flow sub-branches for SP2 for 
H=5 T with the current axis normalized by the critical bias 
current. Splitting of the low-bias Josephson-vortex-flow curve 
into multiple sub-branches indicates the effectiveness of the 
inter-junction capacitive coupling, (b) Josephson vortex-flow 
sub-branches for SP2 at H—4.35 and 4.5 T. The contact re- 
sistance was subtracted numerically in (a) and (b). Inset of 
(a): magnetic field dependence of the capacitive coupling pa- 
rameter, a 



for if=4^5 T [see the inset of Fig. 5(a)], which is con- 
sistent with the assumption of the capacitive coupling 
and reasonably close to the theoretical expectation^ of 
0.1< a <0.4. Therefore, as shown in the upper inset of 
Fig. 4(a), we consider that the first (n—N; the leftmost) 
and last (n=l; the rightmost) sub-branches correspond 
to (i) the triangular and (ii) the rectangular Josephson 
vortex lattices along the c axis, respectively. 

Fig. 5(b) shows the change of Josephson-vortex-flow 
sub-branches for SP2 when magnetic field increases from 
H=4.35 T to 4.5 T. The more changes take place in the 
voltages between the two magnetic fields for branches in 
the higher voltages (branches more to the right). This 
behavior can be compared to the dispersion relation of 
Eq. (1), which states that the lower-index modes are 
more dispersive than the higher-index modes. Thus, the 
right-most branches among multiple Josephson-vortex 
sub-branches in Figs. 3—6 are again considered to repre- 
sent the square vortex lattice [see the inset of Fig. 4(a)] 
corresponding to the n—1 plasma excitation mode. 

The maximum number of vortex-flow sub-branches 
counted for various fields in SP2 is 20, which is two modes 
smaller than the number of quasiparticle branches (or, 
equivalently, than the number of junctions). Recently, 
Ryndyk 34 proposed that the additional dissipation due 
to the charge-imbalance relaxation due to quasiparticle 
injection in a junction can prevent the structural trans- 
formation of Josephson-vortex lattice along the c axis. In 
our sample geometry quasiparticles are injected directly 
from outside through the outermost two junctions, which 
may cause the non-equilibrium charge- imbalance effect in 
the two outermost Cu02 bilayer. The resulting suppres- 
sion of the configurational transformation of Josephson- 
vortex lattice reduces the number of observable resonance 
modes. For inner CuC>2 planes, on the other hand, with 
the pair tunneling for a bias below the critical current of 
the junction, the quasiparticle injection can be neglected. 



B. Temperature dependence 



ence of any defects in the stacked junctions or pancake 
vortices due to any field misalignment, with the converg- 
ing current point corresponding to the depinning current. 
The inductive-capacitive hybrid coupling model indicates 
that the extent of the spread in the low-bias slope of the 
resonant branches is proportional to the strength of the 
capacitive coupling represented by the parameter a as 

V n=[^t] i+2a[i-<jos(r»r/iV)] ' where n=l,2,.. ,N. Here, I p 
is the depinning current of the Josephson-vortex, I c is the 
tunneling critical Josephson current at a given applied 
field. Here V c is the maximum mode voltage correspond- 
ing to a given bias current as represented by the curve 
denoted with an arrow in Fig. 5(a). The gray region 
illustrates the best fit to the sub-branches in the low- 
bias region. The best-fit value of a turns out to be 0.45 
with V c =4.22 mV. The best-fit values of a were found to 
be almost insensitive to magnetic fields as a=0.43^0.45 



Figs. 6(a) and (b) show the temperature dependen- 
cies of I- V characteristics for SP2 in H=5.2 T and for 
SP3 in H=3 T, respectively. The temperature depen- 
dencies of the sub-branches in the Josephson vortex-flow 
region and quasiparticle-tunneling region are much differ- 
ent from each other. Although not clearly illustrated in 
the figures the quasiparticle branches survived up to tem- 
peratures close to T c . The sub-branches in the Josephson 
vortex-flow region, however, disappeared at a character- 
istic temperature T c h, near 18 K and 28 K for SP2 and 
SP3, respectively, leaving a single branch as SP1 in Fig. 
3(b). These characteristic temperatures are almost in- 
sensitive to the fields applied [see insets of Figs. 6(a) and 
(b)]. 

In the vortex flow state for the magnetic-field range of 
Regions 3 and 4 the triangular vortex configuration and 
the multiple-mode configuration resonating with the col- 
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V (mV) V (mV) 

FIG. 6: Temperature dependence of I- V characteristics for (a) 
SP2 in H=5.2 T and (b) for SP3 in H=3 T. The curves are 
shifted vertically for clarity. The insets of (a) and (b) show a 
single Josephson vortex-flow branch above the characteristic 
temperature T c h in two different magnetic fields for SP2 and 
SP3, respectively. 



lective CTP modes compete each other. At low-enough 
temperatures, the triangular vortex configuration wins in 
the low-bias region, while the multiple-mode configura- 
tion predominates in the high-bias region. For T^T c h, 
however, thermal fluctuation and quasiparticle dissipa- 
tions prevent the Josephson vortex lattice from resonat- 
ing with the collective CTP modes. Then the vortices 
stay on the triangular-lattice configuration without a 
structural transformation even in the dynamic stated 
The value of V m seems to be insensitive to temperature 
near T c h, which means that the triangular Josephson vor- 
tex lattice can reach the maximum velocity correspond- 
ing to the fastest mode, n=l, without the structural 
transformation in resonance with the CTP modes. In 
this branch above T c h, therefore, the lattice structure 
may be the modulated triangular lattice but with the 
propagation velocity still locked to that of the fastest 
mode. Thus, the single branch above T c h represents the 
triangular vortex lattice, while the well-defined outermost 
sub-branch below T c h represents the rectangular lattice. 
Accordingly, in Figs. 3(b) and 6, although the rightmost 
sub-branch below T c h appears to directly transform into 
the single branch above T c h, their characteristics corre- 
spond to very different Josephson-vortex configurations. 

As mentioned earlier, the single Josephson vortex- 
flow branch without sub-branches has been observed in 
the usual mesa structured In the usual mesa or the 
S-shaped-stack structure the structural transformation 
from the triangular lattice to the rectangular lattice in 
resonance to the CTP modes is hard to take place be- 
cause of the strong coupling to the near-static triangular 
lattice in the basal layer even at sufficiently low temper- 



atures. In this case, the velocity of the Josephson vortex 
lattice may reach that of the fastest mode without the 
structural transformation. 34 



V. CONCLUSIONS 

The Josephson vortex dynamics was comprehensively 
studied in compactly stacked intrinsic Josephson junc- 
tions over a wide range of magnetic fields and tempera- 
tures. Measurements were made on three stacks of intrin- 
sic Josephson junctions, each sandwiched between two 
(top and bottom) normal-metallic electrodes, thus elim- 
inating the interference from the basal layer (s). The I- 
V characteristics for high external magnetic fields (Re- 
gions 3 and 4 defined in Fig. 1) showed three states: 
(i) the near-static triangular Josephson-vortex state with 
the formation of a single branch, (ii) the dynamic vor- 
tex state with multiple-mode sub-branches, and (iii) the 
McCumber quasiparticle-tunneling state with quasipar- 
ticle sub-branches. In the dynamic vortex state, each 
sub-branch corresponded to the different configuration 
of the Josephson-vortex lattice along the c axis from 
a triangular to a rectangular lattices. 5 With increas- 
ing temperature, the sub-branches disappeared at T c h, 
above which no configurational transformation of the 
Josephson-vortex lattice can be assumed. 

Due to the strong inductive coupling to many junc- 
tions in the basal layer(s), the mesa or the S-shaped- 
stack structure, although containing a finite number of 
junctions in the stack under study, behaves effectively 
as consisting of an infinite number of junctions along 
the c axis. Thus, these structures cannot support the 
formation of the multiple Josephson-vortex-flow modes 
as predicted theoretically for a finite number of junc- 
tions. Although the c-axis boundary condition is not 
strictly imposed in these systems the transverse vortex 
motion is still restricted by the boundary conditions at 
the side edges. This is the main reason why, in the mesa 
or the S-shaped-stack structure, clear magnetoresistance 
oscillations are observed by the interplay between the 
near-static triangular Josephson-vortex configuration 21 
and the boundary potential in the systems, while the 
transverse Josephson-vortex multiple modes are not ob- 
tained. One should be noted that the multiple collective 
Josephson vortex modes can be obtained only in a stack 
containing a finite number of (intrinsic) Josephson junc- 
tions. The triumph of this study is that our samples truly 
satisfied this condition of a finite number of junctions in 
a stack and, as a result, both the near-static Josephson 
vortex state highly coupled to the edge potential and the 
collective transverse Josephson-vortex state were identi- 
fied in each junction. This facilitated the detailed iden- 
tification of the vortex configurations illustrated in Fig. 
1. 
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